Luttinger liquid behavior in tunneling through a quantum dot at zero magnetic field 
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Thermodynamic and transport properties of a two dimensional circular quantum dot are studied 
theoretically at zero magnetic field. In the limit of a large confining potential, where the dot 
spectrum exhibits a shell structure, it is argued that both spectral and transport properties should 
exhibit Luttinger liquid behavior. These predictions are verified by direct numerical diagonalization. 
The experimental implications of such Luttinger liquid characteristics are discussed. 
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Non-Fermi liquid systems - electronic systems whose 
elementary excitations cannot be described by electrons - 
has fascinated physicists due to their unusual properties 
(such as superconductivity and magnetism). Luttinger 
Liquid (LL), describing interacting electrons in one di- 
mension, is one of the most studied models of such a 
non-Fermi liquid system, since it has been solved long 
time ago The non Fermi-liquid characteristics are 
expected to play a major role in transport through a 
LL or tunneling into it giving rise to a power-law 
dependence of the current or the tunneling density of 
states on energy, voltage or temperature. Accordingly, 
there has been many theoretical suggestions of physical 
systems that should exhibit LL behavior and numerous 
experimental attempts to observe this behavior in such 
systems. These attempts have not been very successful, 
due to the sensitivity of LL to disorder, and the difficulty 
of making clean one-dimensional systems. Only recently, 
several decades after the model was originally introduced, 
experimental observation of Luttinger liquid behavior in 
one-dimensional systems such as nanotubes Q and semi- 
conductor systems Q| has been reported. Two dimen- 
sional systems are also expected to exhibit LL behavior 
in strong magnetic fields, due to the formation of one 
dimensional edge-states near the edge of the systems [|j . 
However, tunneling experiments into the edge of a two di- 
mensional electron gas in strong magnetic fields ||^ have 
not been conclusive. 

Quantum dots have been used in the last several years 
to study various aspects of strongly correlated systems, 
most notably the Kondo effect [^jsj. The experimental 
control over the properties of the dot allow a more de- 
tailed investigation of the strongly correlated states than 
that accessible in bulk systems. An experimental study 
of LL behavior in quantum dots has been suggested by 
Kinaret et al. |^ who studied theoretically tunneling into 
a quantum dot in a strong magnetic field. Here, due to 
the finite number of electrons in the dot N , the power-law 
behavior of the tunneling density of states as a function 
of energy, characteristic of LL, will be manifested in the 
power-law dependence of the tunneling probability on N . 



Technological progress allows now a construction of cir- 
cularly symmetric quantum dots ||ict| . The spectrum of 
these dots displays shell structure expected from such 
symmetry. In this letter we propose and demonstrate 
a new type of Luttinger liquid behavior in such circu- 
lar two-dimensional quantum dots at zero magnetic field, 
which is manifested in the spectral and transport prop- 
erties of the dot. The main idea in this work is that 
because of the shell structure, one of the two quantum 
numbers needed to characterize an electronic state in two 
dimensions is frozen (for excitations within the shell), 
and thus the excitations are described by an effective 
one-dimensional system. We explicitly solve the quan- 
tum dot Hamiltonian with long-range electron-electron 
interactions and obtained the low-energy spectrum and 
wave-functions of this many-body system. The spectrum 
and the tunneling amplitude indeed follow the predic- 
tions of Luttinger liquid behavior. 

We describe a quantum dot by a system of J\f inter- 
acting spinless electrons in two dimensions, confined by 
a parabolic potential. 
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with the interaction potential V{r). The noninteracting 
part of the Hamiltonian is trivially diagonalized, and its 
eigenstates and energies are given by 

^„Ar) = Arl^le-'''/2+^^^ iFi{-{n - \e\ - l)/2, \i\ + 1, r^) 

En,e = nhujo, (2) 

with n — 1,2, . . . and i = —{n — 1), — (n — 3) . . . , n — 1, 
and where iFi is the confiuent hypergeometric function. 
r and 9 are, respectively, the radial and azimuthal co- 
ordinates of the two-dimensional vector r, and r is in 
units of yjfi/ mujQ . Thus the spectrum consists of equally 
spaced shells. The energy of the n-th shell is nhcuQ, and 
its degeneracy is n. Within a given shell the states are 
characterized by a single quantum number, e.g. the an- 
gular momentum £. 

The average mean-field Coulomb interactions between 
electrons leads to the usual Coulomb blockade, but does 
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not affect the spectral properties for a given electron 
number A/". For the rest of the paper, while including 
the full Coulomb interaction (the last term in (||)), we 
make the assumption that the inter-shell gap Hloq is sig- 
nificantly larger than the fluctuations in the Coulomb 
interactions between electrons in different levels, which 
are expected to be small Thus the Coulomb inter- 
actions cannot excite electrons from one shell to another. 
Consider now a quantum dot containing Af electrons oc- 
cupying Ush shells. All shells except the last will be fully 
occupied, leaving N = Af — Ushirish — l)/2 electrons in 
the partially occupied shell (POS). Since the Coulomb in- 
teraction do not excite electrons between the shells, then 
the ground state and the lowest excited states will be- 
long to a subspace containing ("j^*") states, each consists 
oi Af — N electrons fully occupying the lower Ugh — 1 
shells and N electrons occupying N of the nsh states in 
the POS. All these states are degenerate in the noninter- 
acting limit. Our aim is to diagonalize the full interacting 
Hamiltonian (|l|) in this subspace. 

In principle, the Coulomb interactions scatter two elec- 
trons from two initial states onto two final states, with 
the conservation of angular momentum. In the absence 
of inter-shell excitations, these terms can be divided into 
three contributions - (a) those where the two initial states 
(and thus the final states) are in the filled shells, (b) those 
where one of the initial states is in the filled shell (and is 
necessarily one of the final states), and (c) those where 
the two states are in the POS. In the relevant subspace 
the first contribution is a constant, and thus is disre- 
garded. The second type of terms, the effective (Hartree) 
potential due to the filled shells , contribute to the di- 
agonal of the Hamiltonian matrix, leading to the removal 
of the energy degeneracy between the different angular 
momenta states (angular momenta ±£ are still degener- 
ate). The resulting dispersion turns out to be very close 
to an inverse parabola. 

Thus, the effective Hamiltonian in the relevant sub- 
space can be written as 
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with 



^ dVrfV|*„.,(r)|2y(r~r')|*„,,/(/)P 

xy(r-r')*„.„,(r)*„,,^,(r') (4) 

where q annihilates a particle with angular momentum 
£ in the POS, and the summations in Tie// are over all 
angular momenta states in the POS. 

This effective Hamiltonian in the POS resembles the 
Hamiltonian describing interacting electrons in one di- 
mension, 
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where pk = '^k+j'^^, Cfe = vpk, with k the momentum 
quantum number and vp the Fermi velocity. Compar- 
ing the effective Hamiltonian (||) to (||) we note that in 
our case the interaction V depends on three quantum 
numbers, and not only on their difference. Thus unlike 
the LL model, the interaction is not merely a product 
of densities. This difference arises because the system 
is not strictly one dimensional, but the electronic wave- 
functions in the Ush — th shell, while having the same 
value of (r^) oc n^/j, do have a nontrivial dependence on 
r. Additionally, the number of the electrons in the POS, 
which are described the the effective Hamiltonian (|^), is 
finite (and could be small). 

Since the one-dimensional Hamiltonian (|^) is exactly 
solvable and is described by LL theory, it is tempting to 
check whether the two-dimensional behavior of a sym- 
metric quantum dot is governed by LL physics. In order 
to investigate the relevance of LL behavior to the spec- 
tral and transport properties of the dot, we diagonal- 
ized numerically the Hamiltonian (^ for several values 
of Ush and N, and for a particular choice of interaction, 
V{r) = log(l/|r|). This choice allows us to evaluate all 
the matrix elements ei and Vj^k,e analytically. To in- 
vestigate the particular role of the interaction term (the 
second term in (^), we multiplied it by a prefactor a 
and carry out our study for several values of a. The 
results shown here are for Ush = 20. For each value of 
total angular momentum L, the Hamiltonian was diago- 
nalized in the subspace of the lowest 6 eigenstates 1 13 , in 



order to determine the lowest energy state with angular 
momentum L. 

For an even number of electrons, the lowest energy 
state is of total angular momentum L ~ 0, the first 
excited states carry angular momentum L — ±2, then 
L = 0, ±4, etc. According to LL theory (for spinless 
electrons), the low-energy spectrum of (H) is given by 
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In order to check the spectrum of the two-dimensional 
quantum dot (^) against the predictions of LL theory, 
we compare in Fig. 1 the evaluated gap, A = {El=2 ~ 
-El=o)/2, to the expectation from LL theory (||). As can 
be clearly seen from the figure, the nontrivial dependence 
of the gap on the interaction strength, a, is very well fit- 
ted by the functional form expected from LL theory, with 
a single fitting parameter Vp=2 — 0.055. (Note that in 
our case the Fermi velocity, vp, is also a-dependent, due 
to the diagonal interaction terms.) Numerically we in- 
deed find that the interaction integrals are about 
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0.05 ± 0.01 for several values of k and £ (independent on 
what branch k and I are). The interaction integrals for 
higher values of p are much smaller, e.g. V4,j^e — 0.01. 
For a = 20, however, we see deviations from the theo- 
retical predictions, indicating perhaps that such strong 
interactions involve states beyond the linear dispersion 
regime. 
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Fig. 1. Comparison of the numerically evaluated gap A in 
the quantum dot spectrum to the prediction of Luttinger 
liquid theory [Eq.(^]. There is a single fitting parameter, 
yp=2 = 0.055 (see text). Inset: Plot of (A^ - Vp)/vF vs. 
the interaction strength. Luttinger liquid theory predicts a 
straight line, of slope 2V2/n. 

LL characteristics are manifested not only in the spec- 
trum, but also in the matrix elements, or the tunneling 
probabilities. The amplitudes of the Coulomb blockade 
peaks, Tjv, in the linear response conductance through 
the quantum dot are proportional, when the tempera- 
ture is smaller than the gap in the excitation spectrum, 
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where ^I^jv is the many-particle ground state of a sys- 
tem with N particles is the POS, and annihilated 
a particle with the corresponding angular momentum 
(namely, the difference in total angular momentum be- 
tween ^'at and "^N+i)- This quantity was studied ex- 
tensively for the case of quantum dot in the fractional 
quantum Hall regime by Kinaret et. al. Using the 
explicit form of the edge-state propagators they showed 
that Tjv oc Ar-(i/'^-i)/2^ where i/ is the filhng factor. This 
dependence was also verified numerically. The fractional 
quantum Hall liquid turned out to be a special case of the 
general theory of tunneling into a LL where it was 
shown that the tunneling probability depends on energy 
in a power-law manner. 



where g characterizes the strength of interactions in LL 
theory (17 = 1 no interactions, g < I repulsive inter- 
actions). The power-law dependence on e is directly 
translated into a power-law dependence of Tn on N, 
Tn oc iV-(/5-i)/2 (the factor 2 comes from the two- 
dimensionahty, N L"^ fc~^ ~ e~^)- 

To check these LL predictions we have evaluated nu- 
merically the matrix elements | < M\dg\'^ n+i > for 
the same set of parameters mentioned above, and for dif- 
ferent values of the interaction strength a. As can be 
seen in Fig. 2, these matrix elements exhibit power-law 
dependence on the number of particles in the POS, in 
agreement with LL theory. Surprisingly the power-law 
behavior is already manifested in tunneling into a quan- 
tum dot with a small number of electrons in the POS. 
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Fig. 2. The numerically evaluated tunneling matrix elements 
Tat, which are proportional to the Coulomb blockade peak 
amplitudes, for several values of the interaction strength a. 
In agreement with Luttinger liquid theory, a power-law de- 
pendence of Tat on A'' is observed, with an exponent that 
increases with the strength of the interactions. 

To further check the relevance of LL theory to the tun- 
neling matrix elements, we employ the functional depen- 
dence of the LL parameter g on the interaction strength, 
which can be explicitly evaluated for a one-dimensional 
model with long-range interactions (only Vq 7^ in (||)). 
For this case ITl , 
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In Fig. 3 we compare the resulting exponent to the pre- 
diction of LL theory with such long range interactions. 
We see an excellent agreement between the numerically 
obtained exponent and the theoretical prediction, with a 
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single fitting parameter Vq ~ 0.33 (with again a deviation 
at very large interaction strength, a = 20). The agree- 
ment between the theoretical prediction and the numer- 
ical calculation also indicates that the logarithmic inter- 
actions used in the numerics mimic the long-range nature 
(V^ = 0,p > 0) of the prediction (^). 
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Fig. 3. The numerically evaluated tunneling matrix expo- 
nent /3, derived from Fig. 2, as a functions of the interac- 
tion strength, a. The dependence upon a exhibits a good 
agreement with the predictions of Luttinger liquid theory, 
Eqs.(|,|). 

To conclude, we have demonstrated, by heuristic argu- 
ments and by numerical calculations, that the excitations 
in a circular quantum dot at zero magnetic field are de- 
scribed by Luttinger liquid theory. These should have 
an observable influence on transport through quantum 
dots. We predict that in a given shell the amplitude 
of the Coulomb blockade peak decreases in a power-law 
fashion. Taking into account that this power-law was 
observed numerically for a quantum dot with a small 
number of electrons in the outer shell (see Fig. 2), this 
prediction can be tested with existing structures. 

An even more intriguing possibility is the relevance of 
this work to transport through a system of interacting 
electrons in the presence of potential disorder, one of the 
most challenging problems of condensed matter physics. 
When the disorder potential is large enough to break the 
electron liquid system into puddles, or effective quantum 
dots, the overall dissipative resistance may be determined 
by transport through such a random quantum dot array 
p6[ . Thus the present calculation suggests that trans- 
port in a two-dimensional disordered interacting system 
may be dominated by LL behavior, similar to transport 
in such a system in a strong magnetic field . 

In this paper the spin of the electrons was ignored. It 
would be quite interesting to investigate the role of spin 



in such devices and explore the possibility of e.g. spin- 
charge separation and its relevance to transport through 
a quantum dot. Another avenue worth exploring is the 
question whether other systems exhibiting shell struc- 
ture, such as nuclei and atoms, can also, under some 
circumstances, exhibit Luttinger liquid behavior. These 
questions will be investigated in a future publication. 
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